TO WHOM IT MAY CONCERN: 



BE IT KNOWN THAT I, DAVID MARK BLUNT, a citizen of 
Australia, residing in South Melbourne 3205, Victoria, in 
the Country of Australia, have invented a new and useful 
improvement in 



SYNCHRONOUS AVERAGING OF EPICYCLIC SUN GEAR VIBRATION 



la 



SYNCHRONOUS AVERAGING OF EPICYCLIC 
SUN GEAR VIBRATION 



This application claims priority from 
provisional application Serial No. 60/442 , 282 , filed 
January 24, 2003 . 

This invention relates generally to the 
processing of vibration data obtained from epicyclic 
gear trains, and more particularly to a method of 
performing synchronous averaging of sun gear vibration 
or vibrations. 

The multiplicity of identical planet gears in 
an epicyclic gearbox, combined with the fact that the 
planet gear axes move in relation to the gearbox casing, 
make traditional vibration analysis techniques 
ineffective in the diagnosis of planet gear faults. The 
inability to effectively diagnose faults in epicyclic 
gearboxes is of particular concern in helicopters, most 
of which use epicyclic reduction gears as the final 
drive to the main rotor. 

As most helicopter transmissions, and many 
aircraft propeller reduction gearboxes, utilize at least 



one epicyclic gear train, and since this forms a non- 
redundant critical part of the drive to the main rotor, 
or propeller it is important to have advanced techniques 
and tools to assess the conditions of these components. 
One such tool is vibration analysis. However, epicyclic 
gear train vibrations are difficult to analyze. Not 
only are there multiple planet gears producing similar 
vibrations, but there are multiple and time-varying 
vibration transmission paths from the gear mesh points 
to any vibration transducers mounted on the gearbox 
housing. 

Previous work resulted in the development of a 
signal processing algorithm that "separated" the 
vibration from individual planet gears in an epicyclic 
gear train. This algorithm proved to be more effective 
at detecting faults in planetary gears than earlier 
planet -separation algorithms, or the ordinary 
synchronous averaging process, but it was not extended 
to the analysis of sun gear vibration at that stage. 
There is need for sun gear averaging method or methods 
that utilize a planet separation algorithm or technique 
to improve the detection of sun gear faults. 



SUMMARY OF THE INVENTION 



It is a major object of the invention to 
provide an improved method of detection of sun gear 
faults. Basically, the technique computes "separated" 
averages of the sun gear vibration "seen" through each 
planet, which are then phase shifted to account for 
their difference in meshing positions and recombined 
(averaged) to produce a modified sun gear average. With 
separation of the vibration from each planet, the 
technique produces a modified sun gear average that 
represents the average meshing behavior of the sun gear 
seen through a single "average" planet. In contrast, an 
ordinary sun gear average computed without planet 
separation represents the average meshing behavior of 
the sun gear seen through all planets simultaneously 
(with the accompanying summation and cancellation of 
certain vibration components) . 

Steps of the new and preferred method include: 

a) detecting sun gear vibrations transmitted 
through each planet gear, 

b) computing separated averages of such 
detected vibrations. 



c) phase shifting said averages to account 
for the differences in gear meshing positions, 

d) and re-combining said phase shifted 
averages to produce a modified average value of the sun 
gear vibration. 

As will be seen, a detection transducer is 
typically provided and operated on the ring gear, and 
sun gear vibration is transmitted to the transducer 
through the individual planet gears, and the expected 
sun gear vibration signal detected by the transducer is 
the sum of the sun gear vibration with each planet 
multiplied by the individual planet pass modulation. 
The measured transducer signal is separated into the 
components from each planet by applying a filter 
consisting of a window function for each planet to the 
measured signal. Sun gear averages are computed from 
each of these separated signals starting at the same 
instant in time. One average is selected as a 
reference, and the remaining averages are phase shifted 
by the angle between its planet and the reference planet 
so that the beginning of each phase shifted average 
starts with the same sun gear tooth in mesh with each 
planet . 



The new technique was applied to vibration 
data collected from sun gear seeded- fault tests in both 
a simple industrial three-planet epicyclic gearbox, and 
a far more complex helicopter gearbox incorporating a 
five-planet epicyclic gear train. The technique 
significantly improved the detection of an isolated 
(single tooth) sun gear fault in both gearboxes. In 
both cases the fault condition index was increased by up 
to 66%. This increase was of more consequence for the 
helicopter gearbox, where the index was improved from 
around the warning level of 3.5, and hence a marginal 
level of fault detection using synchronous averaging, to 
well over the danger level of 4.5 using the new 
technique. It also allowed the fault to be detected at 
lower torques, and improved the ability to use one ring- 
gear accelerometer signal to confirm the detection of 
the fault with another ring-gear accelerometer. 
Accordingly, the new sun gear averaging technique 
significantly improves the condition assessment of 
epicyclic gear trains using vibration analysis. 

These and other objects and advantages of the 
invention, as well as the details of an illustrative 
embodiments, will be more fully understood from the 
following specification and drawings, in which: 



DRAWING DESCRIPTION 



Fig. 1.1 is a view of an epicyclic gear train, 
with vibration detector; 
5 Figs. 1.2 a) and b) are graphs of window 

functions; 

Fig. 2.1 is a view of an epicyclic gear train, 
showing phase shifting; 

Fig. 3.1 is a block diagram showing processing 
10 of detected vibration; 

Fig. 4.1 is an exploded view of a Brevini 

gearbox; 
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Fig. 4.10 is a perspective and schematic view 
of a drive train for a helicopter employing epicyclic 
gearing; 
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Fig. 4.3.1 is a perspective view of a sun gear, 
with one -third of a tooth removed; 

Fig. 4.12 is a perspective view of a 
transmission, showing sensor locations; 
5 Fig. 4.13 (a) and (b) are graphs of window 

functions; 

Fig. 14.14 (a) and (b) are graphs plotting 
residual Kurtosis values as a function of engine torque; 
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INTRODUCTION 



As most helicopter transmissions, and many 
aircraft propeller reduction gearboxes, utilize at least 
one epicyclic gear train, and this forms a non- redundant 
critical part of the drive to the main rotor, or 
propeller, it is important to have advanced techniques 
and tools to assess the condition of these components. 
One such tool is vibration analysis. However, epicyclic 
gear train vibrations are difficult to analyse. Not 
only are there multiple planet gears producing similar 
vibrations, but there are multiple and time-varying 
vibration transmission paths from the gear mesh points 
to any vibration transducer mounted on the gearbox 
housing. 

Previous work at DSTO by Forrester (1, 2 & 3) 
resulted in the development of an algorithm that 
produces "separated" synchronous vibration averages of 
the individual planet gears in an epicyclic gear train. 
This is achieved by applying a window function to be 
vibration signal from a transducer mounted on the ring 
gear. The function weights the signal according to the 
azimuth of the planet gear relative to the transducer 
location. The windowed vibration signal is then 
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synchronously averaged in the usual way. (However, note 
that in order for the separated averages to be uniformly 
weighted, the number of periods averaged must be a whole 
multiple of the number of teeth on the ring gear so that 
there are a whole number of cycles of both the window 
function and the planet gear rotation) . Separation of 
the planet gear averages is achieved by offsetting the 
window function by the angular spacing between the 
planet gears. In effect, this process is attempting to 
use the planet-pass modulation to separate the vibration 
from each planet. Two window functions that have been 
found to perform well are: 



where is the rational frequency of the planet carrier, 
and P is the number of planets (1) . An example of an 
epicyclic gear train with three planets is shown in 
Figure 1.1, and the planet -separation window functions 
for the gear train are shown in Figure 1.2. Note that 
the summation of the individual planet window functions 
is a constant, and therefore the sum of the separated 
averages will be identical to the ordinary {non- 



15 





(2) 
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separated) synchronous average multiplied by that 
constant . 

Forrester's algorithm proved to be more 
effective at detecting faults in planetary gears than 
earlier planet -separation algorithms, or the ordinary 
synchronous averaging process (1, 4, 5 & 6) , but it was 
not extended to the analysis of sun gear vibration at 
that stage. The invention disclosed herein concerns a 
new sun gear averaging technique that utilizes the 
planet -separation algorithm to improve the detection of 
localized sun gear faults. The technique computes 
"separated" averages of the sun gear vibration "seen" 
through each planet, which are then phase shifted to 
account for their difference in meshing positions and 
recombined (averaged) to produce a modified sun gear 
average. Ideally, with complete separation of the 
vibration from each planet, the technique produces a 
modified sun gear average that represents the average 
meshing behavior of the sun gear seen through a single 
"average" planet. In contrast, an ordinary sun gear 
synchronous average computed without planet separation 
represents the average meshing behavior of the sun gear 
seen through all planets simultaneously (with the 
accompanying summation and cancellation of certain 



vibration components) . However, as complete separation 
is impossible in a operational gearbox due to: 

a) the imperfect ability of the window function to 
make full use of the planet -pass modulation to 
separate the planet vibration signals, and 

b) the many vibration transmission paths present 
that are not affected by the planet-pass 
modulation (e.g., through the sun gear shaft, 
bearings, support structure, and housing to the 
ring gear) , 

the new technique produces a modified sun gear average 
that lies somewhere between the ideal case and an 
ordinary sun gear average. 

To test the new technique it was applied to 
vibration data collected from seeded sun gear fault 
tests in both a simple industrial three-planet epicyclic 
gearbox, and a far more complex helicopter gearbox 
incorporating a five-planet epicyclic gear train. 

2. PROCEDURE 
The new procedure or method for computing the 
modified sun gear average contemplates one or more of 
the following: 

1. Use the planet separation algorithm to compute 
the separated sun gear averages. This is almost 
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the same as using the algorithm to compute the 
separated planet gear averages, however the sun 
gear ratio is substituted for the planet gear 
ratio during the averaging process. The same 
separation window function is used as for 
separating the planet gear vibrations (i.e., the 
window function is still based on the planet- 
carrier position) . 

2. Phase shift the separated sun gear averages so 
that the mesh points align. For example, in a 
three planet gear train. Planet 2 must be shifted 
+120°, and Planet 3 must be shifted -12 0° (or 
+240°) to align both with Planet 1, as shown in 
Figure 2.1. Note that the sun gear meshes with 
the planets in the reverse order that they pass 
the transducer on the ring gear. 

3. Combine (average) the aligned averages. 

See for example Fig. 1.1 showing a vibration 
transducer 100 mounted on ring gear 101 of an epicyclic 
transmission 102. The latter also includes three planet 
gears 103 and a sun gear 104. The transducer operates 
as a detector to detect sun gear vibrations transmitted 
through each planet gear, as the sun gear rotates about 
axis 104a. Circuitry at 106 connected at 105 to 100 



12 



operates to compute separated averages of such detected 
vibrations. Circuitry at 108 connected at 107 to 106 
operates to phase shift such averages, as disclosed 
herein, to account for differences in gear meshing 
positions. Circuitry at 110 connected at 109 to 108 
operates to re-combine the phase shifted average values 
to produce a modified average value at 111 of sun gear 
vibrations, displayed or recorded at 112. 

4. THEORY 

The following is adapted from the theoretical 
development found in [1] . 

The meshing frequency, f„, of an epicyclic 
gear train with a fixed ring gear is (refer to Figure 
1.1) . 

where N^, and Ng are the number of teeth on the ring, 

planet and sun gears respectively, and /p and are 

the absolute rational frequencies of the planet carrier, 
planet and sun gears, respectively. The relative 

frequencies (/p + f^) of the planet to the carrier, and 
(/s " /c) of the sun to the carrier are: 
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^ P 

Jm Jc Jc »r 



(4) 



(5) 



Ignoring secondary transmission paths, the sun 



gear vibration is transmitted to the transducer mounted 
on the ring gear through the individual planets. The 
expected sun gear vibration signal measured by the 
transducer will thus be the sum of the sun gear 
vibration with each planet multiplied by the individual 
planet -pass modulations, 



where ap(t) is the amplitude modulation due to planet p^. 

and Vgp{t) is the tooth meshing vibration of the sun 
gear with planet p. Alternatively, this can be 



where 0 is the relative rotation of the sun with respect 
to the planet carrier. 



(6) 



expressed in the angle domain as 




(7) 
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Ignoring any slight variation between the 
planets ; the amplitude modulation function (planet-pass 

modulation), ap{<p), will have the same form of all 
planets, differing only by a phase delay, 27cp/P, 



2np 



V P ) 



(8) 



where a(cp) is the planet-pass modulation function, and 9 
is the planet carrier angle. Since the amplitude 
modulation will repeat with the planet carrier rotation 
period, it can be expressed as a Fourier series, 

where are the Fourier series coefficients. 



The method of extracting the representative 
vibrations of the sun gear "seen" through each planet 
(the "separated" sun gear averages) is to incorporate a 
selective (continuous) time filter into the averaging 
process. The filter proportionally divides the overall 
vibration signal into the estimated contributions from 
each planet. For each separated sun gear average. 
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Zg (9) , the filter window, wO) , is centred at the point 



at which the planet is adjacent to the transducer. 
Assuming that all the vibration that is not synchronous 
with the relative sun gear rotation will tend toward 
zero with the averaging process (7 & 8) , the synchronous 
average taken over N periods of the relative sun gear 
rotation can be expressed as 



P ) 



(10) 



P-\ 



p ){K 



where Vg„(9) is the mean vibration of the sun gear with 
planet k. 

With careful selection of the window 
characteristics and the synchronous averaging 
parameters, the separation can be performed with a 
minimum of "leakage" of vibration from the other planets 
(note that even if the window completely separates the 
vibration of the sun gear with each planet, this will 
not totally eliminate the influence of the sun gear 
meshing with the other planets since there will still be 
transmission paths from those meshes through the sun 
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gear to the planet in question. However, because these 
transmission paths will be more indirect (i.e., the 
vibrations have to travel through the sun gear) , they 
should be less significant.) no loss of vibration data, 
and no distortion of the average. The three conditions 
that meet these criteria are: 

a) the separation window is periodic with 
the planet carrier rotation (i.e., it 
can be expressed as a Fourier series with 

a fundamental frequency of f^) , 

b) the averaging is performed over the 
relative sun rotation period, and the 
number of periods averaged is a whole 
multiple of the number of teeth on the 
ring gear, N^, and 

c) the Fourier series representing the 
separation window has no terms at or 
above the planet -pass frequency {this is 
done to limit aliasing and prevent 
distortion of the average) (i.e.,Np x f^), 
and is a real-valued, even function. 
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Under condition (a) , the separation 
window can be expressed as a Fourier series 



ec 



(11) 



and Equation 10 becomes 




(12) 



Under condition (b) , that the number of 



periods averaged is a whole multiple of the number 
of teeth on the ring gear, jN^ where j is a whole 
number, and since the summation (or integration) of 
any sinusoid over a whole number of cycles will be 
zero, except where the frequency is zero, i.e., 




(13) 



Equation 12 therefore becomes 



Under the condition (c) , that the Fourier 
coefficients = 0 for 1 ^ P, Equation 14 becomes 



hA^) = Z^sA^) ZfV.Aj''-'^ , (15) 



and if the separation window is further limited to 
be a real-valued even function (i.e., Wi = W.^, and 
is real), and if the amplitude modulation is also a 
real-valued even function, this can be reduced to 



,2£ 

p 



(16) 



Note that from Equation 15, the summation 
of the separated averages, z^ pO), can be shovm to 
be equal to the ordinary sun gear average, XgO), 
multiplied by a constant, i.e., 



f=} P'lp-i p-i 



= PiVoA,Y,yM (17) 
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and the separation process is therefore reversible. 

As an illustrative example, if P = 3, the 
first separated sun gear average, using Equation 
16, becomes 



+ V,, (^X^o^ + ^{^'A W,A, cos(f ))) 

+ v,,(^K^, +2{lV,A, cos(f )+ W,A, cos(^))) 



and similarly. 



+ vj0)iW,A,-^2W,A,+2W,A,) 
+ v^A^)iW,A,-W,A,-W,A,) 

v,MW,A,-W,A,-W,A,) ' 
^vJeXW,A,-W,A,-W,A,) 
+ v,,2{0)iF,A, + 2W,A,-^2W,A,) 



Thus, if the amplitude modulation were 
known, and the separation window coefficients were 
set to Wo = l/A^, Wi = l/2Ai, and Wj = I/2A2, the 
separation would be complete. However, this is 
rarely practical in operational gearboxes, and so 
alternative window functions must be used. The two 
separation windows that have been found to perform 
well are those mentioned in the introduction, i.e.. 



Equations 1 and 2, which are repeated here for 
convenience , 



>^poH^(0= T + -cos(2;?r,/) .and w_{/) = - + ^cos(2.'77w/;/). 



Applying the new sun gear averaging 
technique described in §2 to Equation 16, the 
modified sun gear average, Zg„(G), can be expressed 
as 



(18) 



whereby the delay, 27tp/P, aligns the mean sun gear 

vibration with each planet, ,,{0) , so that the 
beginning of each separated average starts with the 
15 same sun gear tooth in mesh with each planet . 

If the sun gear vibrations with each 
planet were identical, they would repeat with a 

period of 27rp/P, i.e., 
v,,(^) = v/^+^l, 

20 ' I PJ (19) 
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where Vg is the meshing vibration of the sun gear 
with a planet. Substituting Equation 19 into 
Equation 18 results in 



z. _ {9) = i7, g + 2f^W,A, cosf l{k - pfdL] 

P-Oi-oL M y P J 



(20) 



and the modified sun gear average would thus 
represent, to within a constant, the average 
vibration of the sun gear with a single planet. 
However, since the sun gear vibration will, in 
practice, not be identical with each planet, the 
modified average will only tend to average-out the 
differences between the meshing behaviour with each 
planet- Nevertheless, because the separated 
averages will be aligned at the same sun gear 
tooth, any localized sun gear defect will always 
appear at the same angular position, and thus be 
reinforced. This should lead to an improved 
ability to detect the fault compared to an ordinary 
sun gear average, where the influence of the defect 
will be more distributed. 
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4 . RESULTS 



4 . 1 Simple Industrial Gearbox 

The new sun gear averaging technique was 
applied to two vibration data sets generated from an 
experimental rig with a small industrial Brevini EMIOIO- 
MN epicyclic gearbox (5) . This gearbox has three planet 
gears with 32 teeth, a sun gear with 28 teeth, and a 
fixed ring gear with 95 teeth (see Figure 4.1). The sun 
gear was driven by a 40 kW electric motor through a belt 
drive, and the planet carrier was loaded by a water 
brake dynamometer through a second belt drive. The 
vibration was measured with an accelerometer mounted 
radially on a small steel block bonded to the outside of 
the ring gear, and a synchronising signal was obtained 
from a photocell on the input shaft. The nominal input 
rotational speed was 19 Hz. The rated power of the 
gearbox at this speed is 18 kW. 

The first data set was measured with the 
gearbox in a no-fault (undamaged) condition. The second 
data set was measured with approximately 0.05 mm ground 
from the loaded face of one sun gear tooth to form a 
narrow flat surface at the pitch line. The gearbox was 



operated at a power level of 12 kW in both cases (i.e., 
67% of the rated power) . 

The sun gear averages were computed over 760 

revolutions of the sun gear with respect to the carrier. 

The results for the no- fault case are presented in 

Figure 4.2 to Figure 4.5, and for the fault case in 

Figure 4 . 6 to Figure 4.9. These figures show the 

following: 

a) Separated Sun Gear Averages ; These are 
the averages obtained by applying the 
planet separation algorithm to the sun 
gear. Effectively, they are the averages 
of the sun gear seen through each 
individual planet gear. Note that the 
separated averages have been normalized 
by scale factors (PW^) obtained from the 
separation windows (e.g., for Wa^„(t) , the 
scale factor is 1.5). 

b) Ordinary Sun Gear Averages ; These are the 
averages computed without the planet 
separation algorithm. 

c) Combined Sun Gear Averages : These are the 
averages obtained by combining the 
separated sun gear averages before 
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aligning the mesh points. Effectively, 
they show that the sum of the separated 
averages is equivalent to the ordinary 
average, thus confirming Equation 17. 
Modified Sun Gear Averages : These are the 
sun gear averages obtained by combining 
the separated sun gear averages after 
aligning the mesh points. 
Residual Signals & Kurtosis Values : The 
residual signal kurtosis is a recognized 
method of detecting a localized gear 
fault. Warning and danger limits for 
this condition index are normally set to 
3.5 and 4.5 respectively. For this 
gearbox, the residual signals were 
obtained by removing the mesh harmonics 
(i.e., multiples of 28 shaft orders) and 
their first and second order sidebands. 
For the fault case, the arrows indicate 
the meshing points of the fault with each 
planet. The kurtosis (or normalized 4^^ 
statistical moment) values for these 
signals are printed on each plot, and 
tabulated in Table 1. Note that, due to 



the dismantling of the gearbox and the 
lack of a planet carrier synchronization 
signal, the numbering of the separated 
sun gear averages may not be consistent 
between the fault and no- fault cases. 
However, this is of no importance to the 
new sun gear averaging technique. 



Table 1: Residual Signal Kurtosis Values 



Separation Window 






Average 


No Fault 


Fault 


No Fault 


Fault 


First Separated Sun Gear Average 


2.S94 


8.540 


2.687 


9.606 


Second Separated Sun Gear Average 


2.673 


5.374 


2.598 


5.918 


Third Separated Sun Gear Average 


2.929 


7.586 


2.868 


8.030 


Ordinary Sun Gear Average 


3.301 


4.870 


3.301 


4.870 


Combined Sun Gear Average 


3.301 


4.870 


3.301 


4.870 


Modified Sun Gear Average 


2.768 


7.860 


2.699 


8.677 



4.2 Complex Helicopter Gearbox 

The new sun gear averaging technique was also 
applied to vibration data from a seeded- fault sun gear 
test in a SH-60B Seahawk main transmission. These data 
were obtained from the United States Naval Air Warfare 
Center (NAWC) through The Technical Cooperation Program 
(TTCP) . The test was conducted in the Helicopter 
Transmission Test Facility at the Naval Air Station at 
Patuxent River, Maryland, in January 2001. This 
facility has the capability to fully load the complete 
transmission including the intermediate and tail 



gearboxes. Power is supplied by the aircraft engines, 
and absorbed by water brakes . 

The SH-60 transmission has an epicyclic stage 
as the final reduction to the main rotor, as shown in 
Figures 4.10. The sun, five planets, and ring gear 
have, respectively, 62, 83, and 228 teeth. The seeded 
fault in this test was the removal of one third of a 
single sun gear tooth, as shown in Figure 4.11. The 
vibrations from the epicyclic stage were measured by the 
port -ring and starboard- ring accelerometers shown in 
Figure 4.12. A tachometer signal from a pick-up on the 
hydraulic drive gear in the starboard accessory gearrbox 
was used as the reference signal for the synchronous 
averaging. The tachometer signal was nominally a sine 
wave with a frequency of 11 kHz at the normal operating 
speed. The ratio of the planet carrier (rotor shaft) 
rotational frequency to the tachometer frequency was: 
3441/8816000. The vibration and tachometer signals were 
sampled at 50 kHz with a bandwidth of 20 kHz. 

The new sun gear averaging technique was 
applied to the vibration data sets detailed in Table 2 . 
The 100% dual-engine power for all versions of the Black 
Hawk and Seahawk helicopters is defined as 2828 shp 
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(2109 kW) , which is equivalent to a torque of 355 lb ft. 
(482 Nm) at the engine drive shaft speed of 20, 900 rpm 
(348.3 Hz). However, the improved durability gearbox 
fitted to the Australian S-70A-9 and S-70B-2 variants, 
and the US Navy SH-60B, actually has a continuous power 
rating of 3400 shp, which is equivalent to a torque of 
427 lb ft. (9) . 

The sun gear averages were computed over 684 
revolutions of the sun gear with respect to the carrier. 
This represented about 60 seconds of data. The 
separation window functions are shown in Figure 4.13. 
The results are presented in Figure 4.16 to Figure 4.27. 
These figures show the following: 

a) Ordinary Sun Gear Averages : These are the 
averages computed without the planet 
separation algorithm. 

b) Modified Sun Gear Averages ; These are the 
averages obtained by combining the 
separated sun gear averages after 
aligning the mesh points. 

c) Residual Signals & Kurtosis Values : The 
residual signals for the Seahawk sun gear 
averages were obtained by removing the 
mesh harmonics (i.e, multiples of 62 
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shaft orders) and their first and second 
order sidebands. A band between 75 and 
110 shaft orders was also removed, as 
shown in Figures 4.15, as this is the 
location of a strong resonance. The 
resonance frequency is slightly higher on 
the port side than the starboard side, 
but the band accommodates this 
difference. The resonance is possibly in 
the gearbox structure or in the 
accelerometer mounting brackets. 
Although the resonance occurs near the 
tail take-off gear mesh frequency, it is 
more likely to be excited by the impacts 
generated by the fault, as these are 
synchronous with the sun gear average 
(and will thus not be averaged- out) . The 
kurtosis values of the residual signals 
are printed on each plot, tabulated in 
Table 3, and plotted in Figure 4.14, 
where the horizontal dotted and dashed 
lines indicate the warning and danger 
limits respectively. 



Table 3: Residual signal kurtosis values. 





Port Ring Accelerometer 


Starboard Ring Acccleromcter 


Torque 
(lb ft) 


Ordinary 


Modified 


Ordinary 


Modified 




bs.JO 


bp^^O 


bs,JO 


425 


3.975 


6.058 (+52%) 


6.186 (+56%) 


4.160 


5.998 (+44%) 


6.900 (+66%) 


350 


3751 


5.451 (+45%) 


6-044 (+61%) 


3.206 


4,863 (+52%) 


5.101 (+59%) 


300 


3.525 


4.288 (+22%) 


3.953 (+12%) 


3.179 


3.966 (+25%) 


4.225 (+33%) 


250 


3.453 


3.682 (+7%) 


3.411 (-1%) 


3.053 


3.840 (+26%) 


3.617 (+18%) 


200 


2.996 


3.929 (+31%) 


4.068 (+36%) 


3.025 


3.805 (+26%) 


3.497 (+16%) 


.165 


2.923 


3.343 (+14%) 


3.090 (+6%) 


3.081 


3.298 (+7%) 


3.148 (+2%) 



5. DISCUSSION 

5.1 Simple Industrial Gearbox 

The results show the following: 

1. The fault is immediately apparent in each of 
the separated sun gear averages, however the 
impact is stronger in the first and third 
separated averages than in the second, as 
shown in Figures 4.6 and Figure 4.7. This 
indicates that there are differences in the 
meshing behaviour of the sun gear with each 
planet. These differences can be attributed 
to physical factors such as slight differences 
between the planets, gear misalignments, etc. 
As the new technique aligns and combines the 
sun/planet mesh points, this effectively 
averages out the differences between the 
planets, and produces an average of the sun 
gear vibration as seen through a single 
"average" planet. 
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The ordinary and combined sun gear averages 
are identical (to within numerical error) . 
Thus, the separation algorithm is reversible, 
and is shown not to distort the vibration 
signal . 

Ordinary synchronous averaging produces a sun 
gear average with three equally spaced fault 
impacts, while the new technique appears to 
produce an average with only one fault impact. 
This is most probably because there is good 
separation of the vibration from each planet, 
and the transducer senses the vibration from 
the mesh point with the planet through which 
the average is taken (i.e., the planet 
immediately adjacent to the transducer) much 
more strongly than the vibration from the more 
distant mesh points with the other planets. 
Note, that three impacts are visible in the 
ordinary average because it is an average 
taken through all the planets simultaneously, 
and the mesh points with each planet are 
distributed instead of being aligned. 



The fault is more easily identified in the 
residual signals of the modified sun gear 
averages than the ordinary averages, as shown 
in Figure 4.8 and Figure 4.9, and Table 1. 
The residual kurtosis jumps from 4.870 to 
7.860 (+61%) using the power separation 
window, and from 4.870 to 8.677 (+78%) using 
the summation separation window. Slightly 
higher residual kurtosis values (8.540 and 
9.606, respectively) are present in the first 
separated sun gear averages for both window 
cases. However, the residual kurtosis values* 
from the modified averages are considered more 
reliable because both window cases also have 
lower residual kurtosis values (5.374 and 
5.918 respectively) for two of the other 
separated averages. Another strong point in 
favour of the new technique is the fact that 
the fault is immediately identifiable in the 
modified average itself, not just in its 
residual signal. Thus, the new technique 
clearly improves the detection of the fault. 
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Complex Helicopter Gearbox 

The results show the following: 

1- The new modified sun gear averaging technique 
significantly improves the detection of the 
sun gear fault using the residual kurtosis 
method. The kurtosis levels are increased by 
up to 66%, as shown in Table 3. This is a 
similar increase to that obtained with the 
simple industrial gearbox. 

2 . The fault is detectable at lower torques with 
the new technique, as shown in Figures 4.14. 
The residual kurtosis levels of the modified 
averages exceed the warning limit of 3.5 for 
all engine torques at or above 200 lb ft., 
with the exception of one measurement (3.411 
and 250 lb ft. for the port-ring accelerometer 
using the bsu„(t) window) . In contrast, the 
residual kurtosis levels of the ordinary 
averages only exceed the warning limit at or 
above 3 00 lb ft. for the port accelerometer, 
and only at 425 lb ft. for the starboard 
accelerometer . 
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The residual kurtosis levels of the modified 
averages surpass the danger level of 4.5 at 
the higher torques of 350 lb ft. and 425 lb 
ft. for both accelerometers, while this level 
is never exceeded with the ordinary averages. 
Thus, the new technique assigns a larger and 
more appropriate condition index level to what 
is a severe fault. 

Apart from the previously noted exception, the 
new technique produces similar results from 
both the port and starboard ring gear 
accelerometers, thereby allowing a more robust 
detection of the fault through cross- 
confirmation of one result with the other. In 
comparison, the results for the ordinary 
averages computed from the starboard ring 
accelerometer do not confirm the existence of 
the fault detected with the port ring 
accelerometer until the torque reaches the 
very high level of 425 lb ft. (120% torque) . 
Thus it is arguable whether this severe fault 
would be detected at all using ordinary 
averaging if secondary confirmation is 
required. 



5. Both separation window functions produce 
similar improvements to the detection of the 
fault, although the summation window, W^^^Ct) , 
is perhaps slightly better at the higher 
torques . 

6. In the modified averages where the fault is 
detectable, the averages exhibit five almost 
equally strong impacts per revolution of the 
sun gear. However, as a result of the 
separation algorithm, it might be expected 
that they should have only one impact, or that 
one should be significantly stronger than the 
rest, as was found in the simple industrial 
gearbox case. The results therefore indicate 
that perhaps there is weaker separation of the 
vibration from each planet gear in this case 
(due to the larger number of transmission 
paths in this gearbox, or the placement of the 
transducers in locations where the planet -pass 
modulation is weak) , or that the severity of 
the fault is overcoming the differences in 
transmission path effects. 

It was noted in §4.2 that the resonance in the 
gear averages between 75 and 110 shaft orders was 
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probably being excited by the impacts generated by the 
fault as it meshed with each of the five planets. 
Supporting evidence for this is shown in Figure 5.1, 
which displays the spectra of the sun gear averages from 
a Seahawk gearbox at the same torque {350 lb ft.) 
without a sun gear fault. While these are from a 
different gearbox, which will necessarily have slightly 
different vibration, it can clearly be seen that the 
amplitudes of the frequencies in the shaded resonance 
band are very low compared to those in Figure 4.15. It 
should also be noted that, while this gearbox did not 
have a sun gear fault, it did have a seeded planet gear 
fault (half a tooth removed form one planet) , which 
would also have produced resonance-exciting impacts as 
it meshed with the sun and ring gears. However, the 
frequency of these impacts would have been non- 
synchronous with the relative sun gear rotation, and 
thus have been attenuated during the averaging process. 
It follows from this, that the amplitude of the 
vibration in the resonance band may therefore be a 
useful indicator of the presence of a sun gear fault. 
However, the seeded- fault in this test is very severe, 
and this method of sun gear fault detection may not be 
so effective with a natural incipient fault. 



6. CONCLUSION 



A new technique has been presented for the 
analysis of epicyclic sun gear vibration using a DSTO- 
developed planet separation algorithm. The technique 
has been shown to significantly improve the detection, 
using the residual signal kurtosis method, of an 
isolated (single tooth) sun gear fault in both a 
relatively simple industrial epicyclic gearbox, and a 
far more complex helicopter gearbox. In both cases the 
kurtosis levels were increased by up to 6S%. This 
increase was of more consequence for the helicopter 
gearbox, whereby the kurtosis was improved from around 
the warning level of 3.5, and hence a marginal level of 
fault detection using ordinary synchronous averaging, to 
well over the danger level of 4.5 using the new 
technique. It also allowed the fault to be detected at 
lower torques, and improved the ability to use one ring- 
gear accelerometer signal to confirm the detection of 
the fault with another ring-gear accelerometer. 

It is concluded that the new sun gear 
averaging technique can significantly improve the 
condition assessment of epicyclic gear trains using 
vibration analysis. 



TABLE 2 



Table 2: Vibration data sets. 
Data File 



IMDSun.l.(X)l 
IMDSun.1.002 
IMD5unl.003 
IMDSun.1.004 
IMDSunl.005 
IMDSun.1.006 



Engine 1 Torque 
(lb ft) 


Engine 2 Torque 
(lb ft) 


Tail Torque 
(lb ft) 


Length 
(seconds) 


425 


425 


800 


360 


350 


350 


800 


360 


300 


300 


800 


360 


250 


250 


800 


360 


200 


200 


800 


360 


165 


165 


1450 


360 
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